Resonant Amplification of Gauge Fields in Expanding Universe 
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We investigate the possibility that gauge fields are amplified in an expanding universe by parametric 
resonance, during the oscillatory regime of a scalar field to which they are coupled. We investigate 
the coupling of gauge fields to a charged scalar field and to an axion. For both couplings, gauge 
fields fluctuations undergo exponential instabilities. We discuss how the presence of other charges 
or currents may counteract the resonance, but we argue that in some cases the resonance will persist 
and that hence this mechanism could have some relevance for the problem of large scale primordial 
magnetic fields. 



PACS numbers: 98.80Cq 

I. INTRODUCTION 

In these last years the concept of parametric resonance 
(henceforth PR) has played a fundamental role in cosmol- 
ogy for the development of the theory of preheating [Q , 
the explosive decay of the inflaton after inflation. There 
have been investigations of several types of fields cou- 
pled to the inflaton: scalars jl| , fermions , gravitational 
waves scalar gravitational perturbations Q], graviti- 

Systems which includes nearly jfj conformal invariant 
fields can be efficiently amplified in an expanding uni- 
verse . This occurs because the Robertson- Walker met- 
ric is conformally related to the flat metric, and therefore 
the equations of motion can be mapped into similar ones 
in Minkowski space-time. In these systems the resonance 
can be efficient without resorting to large couplings ||] (it 
is sufficient to have q ~ 0(1), where q is the amplitude of 
the driving term) . The resonance is not weakened by the 
expansion of the universe and it persists until nonlinear 
effects will shut off it 0. 

It is interesting to investigate the effect of PR on fields 
which are conformally coupled to gravity. For these fields 
particle production due to the expansion of the universe 
is absent B. Therefore PR can be an alternative and 
efficient way to amplify fluctuations of these fields. Ex- 
amples of conformally invariant fields are massless vector 
fields and massless fermions 0. In the case of fermion 
the Pauli blocking is present, but this does not prevent 
the occurring of novel interesting effects Q . For a mass- 
less vector field an exponential Bose enhancement occurs. 
In this Letter we investigate the evolution of gauge field 
coupled to a charged scalar field and to an axion. 



II. SCALAR ELECTRODYNAMICS 

The Lagrangian density for scalar electrodynamics is 
jiven by 



C = --^F^F^ - (D^)*(D^) - VX$*$) (1) 

l07T 



where Z? M = V^ t — ieA^ is the gauge and metric covariant 
derivative, is the gauge potential and F^ is the an- 
tysimmetric field tensor, defined as F^ = V ^A v — V V A^. 
From this Lagrangian the equations of motion are: 



- + 2ieA fi d^ - 



dV 



V V F"» = -Airf + Sire 2 A^\ 2 



e 2 A,,A^<S> = (2) 
(3) 



where j M = ie($<9 M < f>* — < 1>*9 M < I>). We shall consider inho- 
mogeneous linearized fluctuations around homogeneous 
quantities f(t, x) = f(t) + Sf(t, x) (where / is any field) 
in the Robertson- Walker metric ds 2 — —dt 2 + a 2 (t)dx 2 . 
We shall study the evolution of the gauge field in the 
Coulomb gauge V • A = 0. This gauge singles simply out 
the dynamical degrees of freedom for the problem (|l|) 
and its close relation to the gauge invariant formalism 
has been pointed out in p0| . By writing the homoge- 
neous condensate as $(t) = e 1 ®^ p(t) / one has at the 
homogeneous level: 



dV 

3Hp+ — 
op 



-p = 



(4) 



plus an equation of motion for A(£) and the constraint 
Aq = 0/e. We shall consider fluctuations SA^t, x) 
around a background A(t) that is chosen to be zero: 

V 2 

SA T + H8A T 2"<5A T + 47re 2 p 2 <5A T = 47r<5j T (5) 



a 2 
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5Aq = AirSjo - Ane 2 p 2 SA , (6) 



where T denotes the transverse part of a vector and 
is a source term which can be different from zero only 
if one consider quantum or statistical correlation of the 
currents [0 because of the symmetry of the space-time 
background. The constraint equation (^) in Fourier space 
is 



1 



ri-2 



4ne 2 p 2 



SA k = 4:nSjo k 



with Sj ok = e(p8(f>2k - f>84>2k) where Sep = (5<f>i + 
iS<h)/y/2- 

We consider now the homogeneous part of the differ- 
ential equation (^J) in Fourier space: 



SA'^ k + Lu^SAxk = 



(8) 



where u)\ k = k 2 + Ane 2 a 2 p 2 and the prime denotes the 
derivative with respect the conformal time r\ [dr\ — dt/a). 
The general solution y of the inhomogeneous differential 



equation 
method p] 



y = Yh + 4?re 2 



) can be obtained through the Green function 



yi{y')y2{y) - mMy^j[ 
w 



a 2 6} Tk dr}/ 
(9) 



where yt is the homogeneous solution, y\ , yi are the two 
linear independent homogeneous solutions and W is their 
wronskian. 

Equation (||) describes a harmonic oscillator with time 
dependent frequency: during the oscillation of the com- 
plex scalar field Eq. (||) can be reduced to a Mathieu- 
like equation . The solutions to this type of equation 
show an exponential instability cx e Mfc ' ? , for some interval 
of frequencies, called resonance bands. There is a known 
correspondence at the level of equations of motion be- 
tween Eq. (||) and the equation of a massless real scalar 
field \ (rescaled by its conformal weigth) interacting with 
a real scalar field by a term g 2 (f> 2 x 2 - 
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(a X )k = 0, (10) 



where R = 6a" /a 3 is the Ricci curvature. In the case 
of conformal coupling (£ = 1/6) this correspondence is 
exact, while for a quartic potential for cf> (or $) this cor- 
respondence is nearly exact since R ~ 0. The be- 
haviour of parametric resonance depends strongly on the 
time behaviour of the homogeneous scalar field, which on 
turn depends on the form of the potential V($). In the 
following we consider power law potentials as: 



V = A„($*$) 



2« 



(11) 



and we briefly illustrate the different behaviours depend- 
ing on the parameters n and A„, restricting ourselves to 
the case in which the scalar field is the dominant compo- 
nent of energy density in the universe 

For a quadratic potential (n = 1) the driving term 
a 2 p 2 in Eq. (||) decays as rj~ 2 . PR is efficient for 
47re 2 /? 2 3> Ai: the resonance is stochastic and broad, 
and the largest Floquet exponent pk occurs for small 
k fluctuations j^]. In the conformally invariant quar- 
tic case (n = 2) the oscillations of p = ap are given 



by an elliptic cosine and Eq. (m reads as an equation 
in Minkowski space-time (technically speaking it is a 
Lame equation). The resonance structure and the rela- 
tive Floquet exponents pk depend non-monotonically by 
the ratio e 2 /A2 0: long wavelengths k 2 -C A2P0 are res_ 
onant for 1/27T < e 2 /A2 < 3/27T, as well for other values 
JjJ , where po is the initial amplitude for p. For n > 3 the 

driving term a p oscillates and grows m time as 77 2n_1 . 
In this case the resonance for Axk (as well for a\) is very 
efficient @|. 

The case of a simmetry breaking potential 



V = m 2 $*$ + A($*$) 



(12) 



with m 2 < 0, is definitively the closest to the electroweak 
phase transition: $ would play the role of the Higgs field, 
and p rolls down from the zero value (false vacuum) to- 
wards Pmin = iy^m 2 / X (true vacuum). Let us note 
that the most important feature in the case of symmetry 
breaking in this toy model of scalar electrodynamics is 
the Higgs mechanism: the gauge field gets an effective 
mass proportional to the value of the field in the broken 
phase y/m 2 /X. Besides the unpleasant feature of consid- 
ering a massive photon, this fact could completely inhibit 
the resonance in an expanding universe or narrow the 
resonance shifting it to scales proportional to \J m 2 / \. 
However, in the realistic electroweak theory the photon 
does not acquire a mass: the symmetry breaking mech- 
anism gives masses to the W and Z bosons, leaving the 
photon massless. This is an important condition for the 
efficiency of parametric resonance in expanding universe. 
We note that in this case the driving term in Eq. (^) is 
given by the W condensate jl6| . 

Let us discuss for completeness also the behaviour of 
the scalar field fluctuations for the potential ( |l2| ) without 
any constraints on the sign of to 2 . The Fourier compo- 
nents of the real part of the rescaled field fluctuation 
8(j)\ = aS(pi satisfy the following equation: 



'ik 



J ik 



= 



(13) 



where u\ k = k 2 + a 2 m 2 + 3Xa 2 p 2 . The Fourier compo- 
nents of the imaginary part of the rescaled field fluctua- 



tion 



satisfy 



J 2k 



u 2k 



ea 2 (28A ok p' + p5A' ol 



a 
a 

- 3aHpSA ok ) 



(14) 



where ui 2 k = k 2 + to 2 a 2 + Xa 2 p 2 . By using the constraint 
equation (^) and introducing qk = S^k/^T onc obtains 

8ire 2 a 2 
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Ik 



J k2 



2 



a 
a 



-{p' + aHpf 



-~L-2- 



q k = 0- (15) 
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We note that the gauge coupling affects the resonance 
structure of the scalar field: while w\ k is the driving term 
does not determine the resonance bands 



for 



for 84>2k, as would do if the charged scalar field were 
uncoupled to a gauge field p7| . 



Besides the PR effect, one of the two circular polar- 
ization exhibit also a negative effective mass (depending 
on the sign of cf>). This feature could be interesting be- 
fore the stage during which the field coherently oscillates, 
when <fi has a fixed sign. 



III. AXION ELECTRODYNAMICS 



IV. PLASMA EFFECTS 



We consider the effect of parametric resonance for an- 
other type of coupling described by: 



r - — P F^ - -d 

L ~ 167T 2 " 



(16) 



The neutral scalar field can represent an axion in fig] 
(in such a case the vacuum angle 9 is defined to be 9 = 
4>/ /, where / is the Peccei-Quinn simmetry scale) or a 
general pseudo-Goldstone boson in 119]. By considering 
again the scalar field as (f>(t, x) = 4>(t) + 6<f>(t, x) one has 
at the homogeneous level: 



3Hcf) 



dV 







(17) 



and for the transverse degrees of freedom 

■ V 2 -V 
<5A T + HSA T —(5 At + 4irg(f)— x <5A T = . (18) 

In this case it is more convenient to consider the circular 
polarized Fourier component perpendicular to k, A± kl 
whose equation of motion is 



SA'i k + (k 2 ± 4irg</)'k)6 A± k = 



(19) 



In the regime of coherent oscillation of the scalar field 
one has again a Mathieu-like equation where the driving 
term is proportional to the time derivative of the scalar 
field and goes as rf~ x for a power law potential anal- 
ogous to Eq. ([ll]) (always under the assumption that 
the scalar field dominates the energy density of the uni- 
verse Jl4|). The efficiency of the resonance depends on 
the value of the adimensional quantity gcj). The analy- 
sis of Eq. ([l9|), describing the axion decay into photons, 
was probably one of the first application of parametric 
resonance in cosmology j2(J , and it has been recently re- 
analyzed in ]2l|^^ ]. By considering <fi as the QCD axion, 
the smallness of gf ~ 10~ 4 prevents the efficiency of the 
resonance in the case of a potential with temperature- 
dependent mass due to QCD effects [^l). In general the 
resonance is on k ~ Ang(f>u>, where u> is the oscillation 
frequency of </>. However, for V — A<^> 4 /4 and 'iirgf = 1 
we observe a linear growth of A± for k/ u> <C &ngf by 
numerically integrating Eq. ([l7]) and (|l9,). Such growth 
is slower than the quasi-exponential amplification of the 
field fluctuations Scf> on smaller scales M. 



We now comment on effects due to the presence of 
other charged particles Q These effects can be important 
since the universe is considered a good conductor after 
reheating until decoupling. For wavelengths larger than 
the collision length [£3| one could use the Ohm relation 
for the additional current j ac jd = uE. This would add to 
Eqs. flq) and (18) a damping term: 



SA' k ' + ATr(aa)6A' k + uj&T})5A k = . 



(20) 



where we have denoted collectively the transverse compo- 
nent as (5A and with u k their time dependent frequency. 
The time behaviour of conductivity is considered in the 
literature to be a oc T oc 1/a |fl8f and so the term mul- 
tiplying the first derivative in Eq. ( pp| ) is constant: in 
this way the resonance on long wavelengths could be 
washed out if aa is larger than the Floquet exponent. 
For wavelengths smaller than the collision length charge 
separation should be taken into account and a term qual- 
itatively similar to the classical plasma frequency would 
enter as: 



(21) 



where n(j]) and m are respectively the number density 
and the mass of the charged particles. This term would 
play the role of an effective mass, but which would de- 
cay as a(t)~ 3 (which is more rapid than the decay of the 
driving term considered in this paper), leaving open the 
possibility of resonance, in particular for large coupling 
costants. Let us observe that these qualitative estimates 
for conductivity and plasma frequency are obtained by 
thermodynamic considerations. Even in the case of a 
preexisting thermal equilibrium the resonance would am- 
plify the gauge fluctuations. Infact a thermal mass ~ T 2 
for the gauge fiels will not be able to prevent the reso- 
nance and it has the effect of shifting the resonance on 
smaller scales (k 2 must be replaced with k 2 + g e gT 2 a 2 ). 

There is an epoch where the the resonance could pro- 
ceed just as worked out in the vacuum case: in preheating 
after inflation, when it is assumed that the "creation" 



*We are referring to charged fields which are uncoupled from 
the scalar field which oscillates coherently. These fields pro- 
duce currents which are incoherent and are coupled in a 
JfiA^ to the gauge field. 
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of charged particles is contemporaneous to the ampli- 
fication of gauge fluctuations. This amplification will 
last until rescattering and backreaction set in to shut off 
the resonance |3, and terms such as conductivity and 
plasma frequency will appear in this out-of-equilibrium 
process by a self-consistent study of the problem 25 1 (al- 
though their time-dependence will be different from the 
corresponding thermal equilibrium quantities). In this 
case the constraints on the coupling costants e 2 and g 
come only from requirement to protect the potential of 
the scalar field from radiative corrections: this argument 
would give the usual e 2 < 10 -6 for the scalar electro- 
dynamic case. If one requires that the charged inflaton 
gives the right amount of e-folds and CMB fluctuations 
(m ~ 10~ 6 Aip; for a massive inflaton, A ~ I0~ 13 for a 
self-interacting inflaton), then a broad resonance regime 
would be allowed for the gauge field. For the axion elec- 
trodynamic case, the same argument leads to gf ^> I0~ 4 , 
opening a possible new efficient channel for the decay of 
the scalar field <fi into gauge fluctuations. 



V. PRIMORDIAL MAGNETIC FIELD BY 
PARAMETRIC RESONANCE? 

Even if the conformal property of gauge fields in cos- 
mological backgrounds is useful for PR, it is the main 
reason why gauge fluctuations are not amplifield by the 
expansion of the universe Q, in contrast to what hap- 
pens for minimally coupled scalar fields. Therefor the 
conformal invariance of gauge fields is the main issue for 
the explanation of primordial magnetic fields |26[ within 
the inflationary paradigm p8| , [I9| . 

PR could be interesting for the generation of primor- 
dial magnetic fields for two reasons. First, it is intringu- 
ing to note how the exponential growth due to parametric 
resonance is of the same type as the dynamo effect [^6| , 
one of the astrophysical processes postulated to explain 
the observational evidence for galactic magnetic fields. 
As a second point, there are several examples in which 
significant amplification occurs on the maximum causal 
scale allowed by the problem (the coherence scale of the 
field). Since the scale of the magnetic seeds is always 
a crucial issue for a model which aims to explain their 
origin, this feature of PR is very interesting. 

As a straightforward result of section II preheating 
could provide an extra growth in the amplitude o f gaug e 
fluctuations which are produced during inflation fllq , |l9| |: 
for parameters in which long wavelength gauge fluctu- 
ations are amplified, this occurs directly on observable 
scales because of the coherence of the inflaton on the 
particle horizon scale. 

We illustrate how the mechanism could work in a sim- 
ple scalar electrodynamic toy model where the charged 
scalar field is a massless self-interacting inflaton with 
27re 2 = A During inflation, assuming $ in slow rollover, 
a solution for Eq. M) is: 



N 1/2 



HP(-k V ) (22) 



where r\ = — 1/ (Ha) is the conformal time used to model 
a de Sitter era (— oo < r\ < 770 < 0, where 770 is some 
arbitrary time) and the index of the Hankel function is 
so defined 



47re 



2 P 
H 2 



3e 2 M 2 ! 



2A 



(23) 



where the last equality holds during slow-rollover and 
A/ 2 [ is the inverse of the Newton constant G. The solu- 
tion ( |22] ) matches the adiabatic vacuum e~ lkr, /V2k for 
i] — > -00. For v 2 > Eq. fl22]) leads to a slight shift of 
the initial vacuum for the long wavelength fluctuations 
(— krj — ► 0) of the gauge field: 



5A 



Tk 



r(z/) I ' —kr) 



(24) 



By using the definition Bi = ti mn d m A n l a one can relate 
the fluctuation of the magnetic field to the gauge field 
fluctuations 



IB* 



fc 2 |<5A 



Tfcl 



(25) 



Parametric resonance after inflation amplifies the gauge 
fluctuations in Eq. (Q) and |Bfe| 2 in Eq. (25) exponen- 
tially up to a maximum factor I0 12 p4j ]. 

However, this last amplification is not sufficient to ex- 
plain the value required by the observation. Indeed, the 
simple toy model of a self-interacting inflaton is very use- 
ful in order to follow gauge fluctuations at a linear order 
from inflation through preheating, but it misses a super- 
adiabatic amplification during inflation, which was inves- 
tigated by Turner and Widrow |18| . This amplification 
can occur due to an effective negative mass for the gauge 
field during inflation, which can be obtained by breaking 
the conformal invariance through a explicit coupling to 
the curvature, for instance (Ts}| . In string cosmology an 
effective negative mass for the gauge field is generated by 
the dilaton coupling [^7) . 

Instead, as we see from Eq. (pi|), long wavelength 
gauge fluctuations are slightly modified in amplitude and 
spectrum with respect to the initial adiabatic spectrum 
by the coupling with inflaton during slow-rollover. In- 
deed the positive mass generated by the coupling of the 
gauge field to the charged scalar field lead to a slight 
suppression. When slow-rollover is ended and v 2 < the 
gauge modes are suppressed. This phase is just a short 
transient, after which the resonant amplification starts 
immediately (roughly when p first cross zero). The main 
point is the decay of B 2 as 1/ a 4 once a gauge fluctuation 
has left the Hubble radius during inflation, leading to a 
suppression which the following preheating phase cannot 
fill up || . 

Now we quantify our previous statements and we es- 
timate the order of magnitude of the amplitude for the 
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magnetic field at decoupling. We shall consider the effect 
of preheating as a k-independent amplification factor tt 2 
which multiplies the spectrum at the end of inflation (24). 
We also neglect the finite time duration of the preheating 
phase and we identify the end of inflation with the begin- 
ning of the usual radiation dominated Friedmann era. At 
the beginning of the radiation dominated phase (denoted 
in the following by 77*) the spectrum of the magnetic field 
is: 



fc 3 |B fe (7f)| 2 = n s 



fc 5 



a 4 (77*) 



\SA Tk ( V n\' 



(26) 



where we consider SAxkii]*) — l/v2fc, which has been 
obtained from Eq. d24) by setting v = 1/2. By using 
the the relation a 4 B const, for the magnetic field 
evolution in the radiation dominated era, the spectrum 
at the decoupling time //dec is: 



fc 3 |B fc (r7 DEC )| 2 = fc 3 |B fe (77 



rp4 

*\|2 1 DEC 



rp4 

REH 



(27) 



since the temperature T scales as the inverse of the scale 
factor a. By using Eqs. ( p6| ) and ( p7| ) we estimate at 
time of decoupling a magnetic field of size 



B|dec(*0 = (fc 3 |B fc (77 DEC 



I2U/2 



O T 2 
" (fckpc) 2 — ^ 
y/2 kpc 2 T 2 



NOW 



10~ 42 Gauss 



(28) 



where it has been assumed a(ryNow) = 1, a comoving 
wavenumber k = 0.1 kpc -1 = 0.64 x 10~ 36 GeV, fl = 10 6 , 
TdecAnow - 10 3 , 1 Gauss = 3 x 10- 20 GeV 2 The 
estimate ( p8|) is lower than the magnetic seed field re- 
quired by observation, which is 10 _34 G for a flat, low 
density universe PQ] . We note that B is proportional to 
fl, but not to Treh or to the Hubble constant during 
inflation, as observed in pq |. This is due to the the ap- 
proximation v = 1/2 which restablishes the conformal 
invariant spectrum in the estimate (Pq). 



In the simplest toy model with a massless charged in- 
flaton, a linear calculation based on preheating does not 
provide directly an amplitude in agreement with obser- 
vations because the magnetic field fluctuations are red- 
shifted during inflation. However, this latter effect may 
be circumvented in many ways. In more realistic infla- 
tionary models motivated by particle physics the gauge 
field can be sustained during inflation and then amplified 
during preheating. A sustain during inflation could be 
also realized by breaking conformal invariance with ex- 
plicit coupling between the gauge field and the Riemann 
curvature in the lagrangian [lq ] . Conformal invariance is 
also effectively broken by the inhomogeneities genarated 
during inflation (a second order effect in perturbations). 
Moreover, PR would also provide an additional growth 
in the generation mechanism based on the statistical cor- 
relation of the currents studied in |p"l|| . 

A final amplitude and spectrum of gauge field fluctu- 
ations generated by PR is however beyond the scope of 
the present paper. Infact, its calculation would include 
the analysis of nonlinear effects and a simple estimate 
based on matching techniques is inadequate. In models 
with interacting scalar fields, nonlinear effects as rescat- 
tering pi| play an important role in the generation of 
the final power spectra. Indeed, such power spectra, as 
|<5Afc| 2 ~ k~ a with a > 1, would be very interesting 
for the problem of large scale magnetic fields, because 
steeper than (|24|) in the infrared region. 
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VI. CONCLUSIONS 



In this paper we have studied the possibility of an en- 
hancement of gauge fields due to parametric resonance. 
For both the couplings to a scalar charged field and to 
an axion, the resonance can counteract the redshift due 
to the expansion of the universe and exponential grow- 
ing modes appear. The resonance is generically on small 
scales in the axion case, while it can occur for scales up to 
the coherence scale of the scalar field in the scalar electro- 
dynamic case. Therefore the latter scalar field coupling 
is more interesting for the problem of cosmological mag- 
netic fields, in which the large scale observed seems a 
crucial issue. 
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